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ON THE COMPLETELY INTEGRABLE CALOGERO TYPE 
DISCRETIZATIONS OF NONLINEAR LAX INTEGRABLE DYNAMICAL 
SYSTEMS AND THE RELATED MARKOV TYPE CO-ADJOINT ORBITS 

ANATOLIJ K. PRYKARPATSKI 


Abstract. The Calogero type matrix discretization scheme is applied to constructing the Lax 
type integrable discretizations of one wide enough class of nonlinear integrable dynamical systems 
on functional manifolds. Their Lie-algebraic structure and complete integrability related with 
co-adjoint orbits on the Markov co-algebras is discussed. It is shown that a set of conservation 
laws and the associated Poisson structure ensue as a byproduct of the approach devised. Based 
on the Lie algebras quasi-representation property the limiting procedure of finding the nonlinear 
dynamical systems on the corresponding functional spaces is demonstrated. 


1. Introduction: the discretization and related Markov algebra splitting 

With a fairly generous definition, a one-dimensional real-valued discrete nonlinear dynamical 
system on a manifold M C ^(^R" 1 ) for some finite m G Z + is any evolution equation that 
can be written down as 

(1.1) du/dt = K[u \, 

where t G R is the evolution parameter, u G M and K : M — > T(M ) is some smooth enough 
vector field mm on the manifold M. Very often such equations m can be naturally obtained as 
the standard discretization IIHH of a given smooth nonlinear differential dynamical system 

(1.2) du/dt = /C[u] 

on a functional submanifold A4 C /^(R; R m ), generated by a smooth vector field /C : A4 —> T(M). 
Namely, there exist such mesh points Xj ^ Xi G R for i / j £ Z, that the corresponding vector 
{u(a;j) G R m : j G Z} = u G M and the suitable discretization of (11.21) coincides with (11.11) . 

Other approach to the discretization of la is based on the Calogero type mm scheme of 
constructing finite-dimensional quasi-representations of the infinite-dimensional Heisenberg-Weyl 
algebra of operators h :={f, D x , 1 : x G R}, where D x := d/dx , in some functional submanifold Mq 
C C°°([a, 6]; R) P\L 2 ([a, 6]; R) of differentiable functions, as owing to the well known von-Neumann 
theorem nzii3si, there exists no exact representation of f) in a finite-dimensional functional subspace 
M$ C Mq for any N G Z+. For example, any smooth scalar function / G Afo on an interval 
[a, 6] Cl can be interpolated [26], ;7j in a polynomial form as follows: 

N 

(1.3) f(x) -A f N {x) := 

1=1 

N 

e j( x ) ~ n o -Xi),pj~ n 

i=l ,N,i^j i=l,N ,i^j 

and its derivative, respectively, as 
(1-4) D x f(x) ->■ D x f N (x) = 

i,j=l,N 
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where /n{x) £ Mq subject to the polynomial basis {ej(x) £ Mq : j = 1,.ZV}. Then the well 
known Calogero type quasi-representation m m m m hb of the Heisenberg-Weyl algebra f) is 
obtained as 

(1.5) End{ Mo) diag{xi, x 2 , ..., aqv} G End K), 

End(A4 q) 3 1 —> I := diag{ 1,1,1} £ End l 2 (Zjy; R) 

End(Mo) 3 D x ->• Z := {Z t j := (a;* - 2 Cj) _1 ,i j = 1,1V; 

N 

Z H := ^ (27 — Xj ) _1 : i = 1, IV} £ End 1 2 (Zn; R), 

where interpolating mesh points Xi ^ Xj £ R, 1 ^ j = 1, IV, are chosen to be different and satisfying 
in a suitably defined finite-dimensional Hilbert space I2 (Zjv;M) the strong as N — >• 00 limiting 
canonical Lie algebra relationship 

(1.6) lim ([Z, X\ — I) = 0. 

N—yoo 

The matrix quasi-representations m make it possible to construct easily a naive matrix 
discretization of the nonlinear dynamical system 03 as follows: 

(1.7) dU {m) /dt = lC{U {m \ [Z( m \ C/ (m) ], [Z^ m \ [Z^ m \ t/ (m) ]], ..\Z^ m \ ..., t/ (m) ]]), 

(p-times) 

where the matrices 

(1.8) t/ (m) := u(X (m) ) = diag(u(x 1 ), ti(s 2 ),..., u(®jv)) £ End l 2 (Z w; R)®™, 

Z (m) £ End l 2 (Zjv;M)® m 

(m-times) 

belong to the tensor product matrix space 

(1.9) Endl 2 (ZN', R)® m := End^^bv; R) ® Endl 2 (ZN', R) < 8 >...<S> End^Zjv; R). 

(m-times) 

When deriving the matrix equation & we took into account that 
/C[u]:= K,(u,D x Vl,D^.u, ...Dv.xl) for some fixed p £ Z+, and for arbitrary operator mapping 
<Pn(%) ■ -Mq Mq we made use of the Calogero type quasi-representation property: 

(1.10) End(Mo) 3 (D^ n )(x) -3 {Z , [Z, [Z,..., [Z,g>(X)]}}...} £ End h(Z N ;^), 

( n-times) 

which holds for arbitrary operator derivatives (£)"y>jv)(x) : ~3 -Mq , n G Z+. 

Now we will take into account the observation m that the quasi-representations m belong, 
respectively, to the Markov direct sum splitting of the general Lie algebra gl(N;R) := g = 
M(g) © E(g) : 

(1.11) 7, A' £ E(g), Z£M(fl), 
where, by definition, the linear subspaces 

(1.12) M(g) := (M(A) £ g : M(H) = A - diag(eA)} 

E(g) := {E(H) £ g : E(A) = diag(eA), e := (1,1,.., 1) £ h(Z N -,R)*, 

are Lie subalgebras of the Lie algebra g. Introduce now, by definition, projections Pm( 8 j : g = 
M(g) C g, Pe( 0 ) : g = E(g) C g. Then within the standard R-matrix approach [3] 0 H2s 23 ’25i 
the expression 

(1.13) [X,Y] n := [R X,Y] + [X,RY], R := l/2(P Mffl ) - P E(fl )), 

for arbitrary X, Y £ g defines on g a new Lie commutator structure, generating on the space 77(g) 
the deformed Lie-Poisson bracket 

(1.14) { 7 , g}n(a) :=< a, [Vy(a), V?/(a)] R >=< a, [P M ( 0 ) Vy(a), P M ( 0 ) Vr/(a)] > 
for 7,7 £ 77(g) and any a £ g*, generalizing the classical Lie-Poisson bracket 

(1.15) { 7 , T]}(a) :=< a, [V 7 (a), Vr/(a)] >=< a, [V 7 (a), V?/(a)] > 
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on g. Here the bi-linear trace-functional on the g 

(1.16) < A, Y >:= tr (AY) 

for A, Y 6 g is nondegenerate and Ad- invariant. Taking into account that with respect to this trace- 
functional (11.161) the Lie algebra g ~ g*, the Poisson bracket (II. 141) generates for any Hamiltonian 
function H £ (Z?(g) the following dynamical system on arbitrary a £ g : 

(1.17) da/dt = P E ( 8) J-[P M (0)V£r(a),a], 

where we took into account that projections Pm( 8 ) — f > E(g)- L and P E (b) — PM(j) 1 ■ This construction 
becomes more simpler in the case when the Hamiltonian function H £ 1(g) is taken to be a Casimir 
one with respect to the classical Lie-Poisson bracket (11.151) , satisfying the condition 

(1.18) [a,VH(a)]= 0 
for any a £ g. 

Amongst nonlinear differential dynamical systems (11.21) there exist a wide class of nonlinear 
evolution equations which are Lax type [EJ [Ill IH CE2! integrable and whose discretizations are 
often very important for their numerical analysis and diverse applications. Yet in general, the 
presented above directly discretized matrix dynamical system m does not a priori inherits the 
Lax type integrability of m- Thus, a natural question arises: how to construct a priori Lax type 
integrable matrix discretization of a given Lax type integrable nonlinear dynamical system im? 

As the Lax type representations of the presumably integrable dynamical systems (11.21) depend 
on an arbitrary spectral parameter A £ C, it motivates us to study their corresponding matrix Lax 
type quasi-representations, also depending on the spectral parameter A £ C, as well as depending 
on the basis matrix representation operators m, belonging, respectively, to the Markov direct 
sum splitting of the general Lie algebra g/(A; R) := g = M(g) ©E(g). This can be done effectively 
by means of introducing the notion of the metrized loop mmmm algebra g := g<g>C[[A, A *]] 
and the related Lax type integrable Poisson flows on it. 

Below I present a solution to this posed above question in the case of a special class of nonlinear 
Lax type integrable dynamical systems on functional manifolds making use of the Calogero type 
discretization scheme and the analysis of the Markov type co-adjoint orbits by means of the related 
Lie-algebraic techniques. 

2. The Lie algebraic setting and the Calogero type linear matrix spectral 

PROBLEMS 

There is introduced a metrized loop algebra g := g®C[[A, A -1 ]], generated by the Lie 
algebra g : = gl(N; R) and the related Laurent series 

(2.1) g : ={X(X) = Y X i X '■ X 3 ^0, Z3j'«oo,AeC}. 

j<^oo 

It is endowed with the standard matrix Lie bracket 

(2-2) [A(A),Y(A)] := Y AS ( E 

s<oo j-\-k=s 

defined for any A (A), Y(A) £ g and the simplest Ad-invariant nondegenerate scalar product 

(2.3) < A (A), Y( A) >_i= Tr(A(A)Y(A)) := res tr(X(X)Y (A)), 
satisfying for all A(A), Y(A) and Z( A) £ g the condition 

(2.4) < A(A), [Y(A), Z(X)] >_i=< [A(A), Y(X)],Z(X) >_r . 

Consider now the introduced above Markov splitting (11.121) of the Lie algebra g : 

(2.5) g := M(g) © E(g), 
where for any A £ g 

(2.6) M(A) := A - diag{e A), E(A) := diag(e A), e := (1,1,..., 1) e 1 2 {^n)*, 

and the components M(g)c 0 and E(g) C g are suitable matrix Lie subalgebras, that is 
[M(g),M(g)] C M(g) and [E(g),E(g)] = 06 E(g). 
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The following observation is crucial for our next analysis: the loop algebra g inherits the Markov 
type splitting (12.51) into two Lie subalgebras: 

(2-7) 0 = 0 +© 0 -, 

where the projection P + g : =g + , 

(2.8) g+ := {X(X) = Y XjX 3 : X 0 £ M(g), Xj £ g, N 9 j « oo, A G C} 

iez + 

and the projection P_g : =g_, 

(2.9) g_ := {X(X) = Y YjX j : Y 0 £ E(g), Yj £ g, j£ Z_\{0} , A e C} 

satisfy the Lie commutator relationships [g+,g+] C g+ and [g_,g_] C g_. It is easy to calculate 
the adjoint spaces 0 + and gl to Lie subalgebras (12.8p and (12.91) : 

(2.10) r+ - 0- = {^(A) = ^ YjX 3 : Y 0 £ E(g)\ Yj £ g, j£ Z_\{0} , A G C}, 

iez_ 

0 _ - 3+ = {X(X) = Y, XjX ■ x o £ M(g)\ Xj eg, N 9 j « oo, A e C}, 

iez + 

took into account that 

M(g)* ~ E(g) x = {Y £ g : diag(Y) = 0}, 

E(g)* — M(g) x = {X £ g : X = q 0 e, 

e:=(l,l,...,l) £ hC^N', R)*, q £ h{^N\ R)}- 

The splitting m makes it possible to define a related classical R-structure on the loop algebra 
g : for any X (A), Y( A) eg the commutator 

(2.12) [X(A),Y(A)] r := ([RX(A),Y(A)] + [X(A),RY(A)]), 

satisfies the Lie algebra commutator property, where the linear space homomorphism R : g —> g is 
defined for an arbitrary X(A) £ g as 

(2.13) RX(A) := 1/2(P+X(A) - P_A'(A)). 

The following important classical (see, for instance, [H 151123|i theorem holds. 

Theorem 2.1. (Adler-Kostant-Souriau) Let smooth functionals 7,77 : g* —>R be Casimir ones 
subject to the Lie bracket \2.‘2\) . that is 

(2.14) [V 7 (/(A)), l( A)] = 0 = [Vt 7 (Z( A)), l{ A)]. 
for any l( A) £ g*. Then their modified Lie-Poisson bracket 

(2.15) { 7 , 77 } :=< l( A), [V 7 (/(A)), Vt 7 (Z(A))]r >_i 
vanishes: { 7 , 77 } = 0 on the whole space g*. 

Based on the splitting m one can easily to calculate the actions of adjoint operators P+ : 
0 —> gl C g_ and PL : 0 —> g 1 ). C g+. Namely, owing to the identification g ~ g* one finds that 
the following equalities 

(2.16) p; = p r : 0 -A gl C g_ 
and 

(2.17) 
hold. 


where we 
( 2 . 11 ) 


PI = Pj; : g g\ C 0+ 
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Theorem [2711 above and the equalities (12.161) and (12.171) make it possible to construct a wide 
class of Liouvile integrable dynamical systems mm on the Markov matrix subspace E(g), if to 
reduce the related Hamiltonian vector field 

(2.18) ja{ A) := { H ,a(A} = [P+Vff (a(A)),a(A)], 

on the element a(A) £ g~ g*, generated by a specially chosen smooth Casimir functional H : 
0 —> R. The latter is considered with respect to the standard Lie-Poisson bracket 

(2.19) {7 Mlu :=<Z(A),[V 7 (Z(A)),V7 ? (/(A))] >_, 

for any smooth functionals 7,77 G 2?(g). 

Consider the following smooth functionals 7 „ : g* — > R, n, k G Z+, where 

(2.20) 7 £ fe) := 1 /{n + 1) < (a(A)A” l “ (A)l ) rl+1 , A fc+|a(A)l >_ x = 

= l/(n + l)Tr(A fe+|Q(A) l (a(A)A- |Q(A) l) n ). 


They are, evidently, Casimir ones for the Poisson bracket (12.191) . whose gradients equal 
(2.21) V 7 i fe) (a( A)) = (a(A)A- | “ (A)l )”A fc . 

Taking, for example, n = 2 and k = 4, the corresponding value of the element P+Vy^ (a(A)) G §+ 
at the element a(A) = A 3 1 + A 2 U + AH + Z then equals 


( 2 . 22 ) 


P+Vt4 4 ) (A 3 / + A 2 U + A V + Z) = A 4 / + 2A 3 U+ 
+\ 2 {2V + U 2 ) + A(2 Z + UV + VU) + M {ZU + UZ) . 


From the commutator relationship (12.181) one obtains that the following matrix equation 


(2.23) 


f t (\ 2 U + A V + Z) = [A 3 / + 3A 2 C/2 + A(3Z/2 - 3U 2 /8)~ 
—C7 3 /16 — 3{ZU + UZ)/8 ,X 2 U + XV + Z] 


holds for any parameter A € C. The trace-functionals H m := tr(\ 3 I + \ 2 U +XV + Z) m / 3 , m G Z+, 
are all nontrivial and involutive with respect to the Poisson bracket (12.151) conservation laws of 
the Riemann type discrete matrix dynamical system 


(2.24) 


dU/dt = -2[Z,V], 


dV/dt = -[Z, VU + UZ] + [UZ + ZU, V], 

easily following from (12.231) . Moreover, the discrete matrix dynamical system (12.241) . as follows 
from is a completely integrable discrete approximation of the corresponding limiting as N —> 00 
Riemann type equations in partial derivatives 

(2.25) du/dt= —2v x , dv/dt = uv x — vu x . 

Thus, concerning the limiting dynamical system (12.251) one can formulate the following proposi¬ 
tion. 

Proposition 2.2. The dynamical system \2.25\) allows for any N G Z + the completely integrable 
matrix discretization 1(2.241 ), which allows the Lax type representation 12.231) . 

As a simple consequence of Proposition 12.21 the obtained above limiting dynamical system 
(12.251) is also Lax type integrable. More applications of the devised above Calogero type discrete 
approximation approach to different Lax type integrable nonlinear dynamical systems are under 
preparation. 

























6 


ANATOLIJ K. PRYKARPATSKI 


3. Conclusion 

An observation, that the Calogero type discretization of the Heisenberg-Weyl algebra is related 
to the Markov type splitting of the general Lie algebra gl(N;M),N G Z+, proved to be both 
interesting and useful for analytical constructing completely integrable discretizations of the Lax 
type integrable nonlinear dynamical systems on functional manifolds. They are represented as co- 
adjoint orbits on the Markov type co-algebras and analyzed by means of the modern Lie-algebraic 
techniques. It is shown that a set of conservation laws and the associated Poisson structure ensue 
as a byproduct of the approach devised. Based on the Lie algebras quasi-representation property 
the limiting procedure of finding the nonlinear dynamical systems on the corresponding functional 
spaces is demonstrated. 
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